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SUMMARY 

A method for adaptive identification of reduced-order models for continuous 

stable SISO and MIMO plants is presented. The method recursively finds a model 

whose transfer function (matrix) matches that of the plant on a set of frequencies 

chosen by the designer. The algorithm utilizes the Moving Discrete Fourier Transform 

(MDFT) to continuously monitor the frequency- domain profile of the system input and 

output signals. The MDFT is an efficient method of monitoring discrete points in the 

frequency domain of an evolving function of time. The model parameters are estimated 

from MDFT data using standard recursive parameter estimation techniques. The 

algorithm has been shown in simulations to be quite robust to additive noise in the 

inputs and outputs. A significant advantage of the method is that it enables a type of 

on-line model validation. This is accomplished by simultaneously identifying a number 

of models and comparing each with the plant in the frequency domain. Simulations of 

the method applied to an Sth-order SISO plant and a 10-state 2-input 2-output plant 

are presented. An example of on-line model validation applied to the SISO plant is also 
presented. 
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1. introduction 


This paper presents a new method of robust adaptive identification for stable 
continuous SISO and MIMO systems via recursive rational interpolation. Interpo at ion 
unlike many identification schemes, is frequency-domain based. It attempts to in a 
model frequently a reduced-order model (ROM), whose frequency response (matrix) 
matches that of the plant at certain frequencies. In this case, the model is said to 
interpolate the plant data at these frequencies. The goal, of course, is to in a mo e 
whose transfer matrix closely matches that of the plant at all frequencies. 

The interpolation approach has been studied by several researchers (e.g. W'W)- 
The method presented in this paper is most similar to |4]. In that reference, the 
problem is cast in the form of a system of simultaneous linear algebraic «>“ at,ons “ 
which the unknowns are the plant parameters. We improve upon the method of [4], 
however, in that the present method is recursive and therefore much more robust to 

noise. 


The algorithm presented in this paper monitors a number of points in the 
frequency profiles of the plant inputs and outputs via a recursive version of the Discrete 
Fourier Transform, which in this paper we call the Moving Discrete 
or MDFT. This transform has been discussed in [4] in conjunction wit m 

reduction The transform was also used in another application in (5], where it was 
called a freguency sampling filter. The MDFT has the advantage that it is recursive 
and therefore gives the evolving DFT at every time step. It has the further advantage 
that it is not necessary to calculate the DFT at all N frequencies (for an N - point DFT), 
as is necessary for standard DFTs and FFTs, but only at those frequences desired, 
which are usually far fewer than N in number. This makes the algorithm numenca y 

efficient. 


The paper is organized as follows. Section 2 briefly discusses the MDFT 
algorithm. Section 3 describes the SISO identification algorithm and gives a numerical 
example in which a 2nd-order model is derived for an 8th-order plant. Section 4 
describes a method of on-line model validation made possible by the algorithm and 
gives a numerical example. Section 5 describes the MIMO extension of the algorithm 
and gives a numerical example in which a 6-state model is derived for a 10-state, 2- 
input, 2-output plant. This ROM is compared to several others for the same plant. 

Conclusions are given in Section 6. 
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2. THE MOVING DISCRETE FOURIER TRANSFORM (MDFT) 

Let {a:,}, * = 0, 1, 2, N -1 be a set of N consecutive samples of a continuous 
signal x(t) sampled every At seconds. The well-known Discrete Fourier Transform 
(DFT) of {x,} is the set of complex numbers {X k }, k = 0,1, 2, N-l defined by ([6]): 

** = £ x i Wk ' - * = 0, 1, 2, N-l (2.1) 


where W = exp(- j2tt/ N). Since there are N points in the calculation, this is called an 
iV-point DFT. Under certain circumstances, the complex number X k can be considered 
as a frequency component of x(t) at the frequency f k = k/NAt Hz. We will call the 
frequencies f k DFT frequencies. 


To make the process recursive, assume that a new sample of x becomes available, 
x N . Our objective is to calculate a new DFT of x based on the latest N samples of x. 
That is, we will create a “moving” DFT by discarding the oldest sample, in this case x 0 , 
and calculating the DFT of the new sequence {x,}, t = l, 2, N. In general, we will be 
interested in calculating the DFT of the (n + l)th sequence of a:, i.e. {x„_^ +2 , ..., z n+1 }, 
recursively in terms of the DFT of the nth sequence of x, i.e. {x n _ iV+1 , ..., x n }. If A**"} 
is the *th frequency component of the DFT of the nth sequence of x, then it" can be 
shown ( see M) that -*fc {n+,} is expressed recursively in terms of X k {n} as follows: 


^ {n+1} = — x n _ JV+1 + x n+1 ] W~ k (2 . 2 ; 

Equation (2.2) gives the algorithm for the 1-step-ahead MDFT. Note that the 
algorithm is very efficient because it involves only the addition of a real number to 
X k { } and multiplication by a complex constant W~ k . Referring to (2.2), note that 
+ } depends only on A* {n} and not on any other frequency component. This 
means that it is only necessary to calculate the MDFT at those frequencies which are 
desired, not all N frequencies , as is necessary in standard FFTs. This can result in 
significant savings in calculations since usually derivation of a reduced-order model (see 

Section 3) requires the DFT at a relatively small number of frequencies, whereas N is 
usually large. 
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3. RECURSIVE SISO REDUCED-ORDER MODELING 


Assume that a fcth-order model for a given stable plant is desired, 
model can be described by a fcth-order differential equation of the form: 

Then, the 

a{D)y(t) = p(D)u{t ) 

(3.1) 

where a{D) is a /tth-order monic Hurwitz polynomial in the differential operator 
D = d/dt, and P{D) is a polynomial of order k-1 or less. That is, 

a(D) = D k +'£a i D i 

i= 0 

(3.2) 

w>) = 

»=0 

(3.3) 

Equation (3.1) can be rewritten in the form: 


D k y = -a fc _! D k ~ l y- . . ■ ~a 0 y + /? fc _ x I>* -1 u+ • - • + Po u 

(3.4) 

which can be rewritten as 


Y(t) = 4’ T Z(t) 

(3.5) 

where 


Y(t) = D k y 

(3.6) 

X /> T = [— • i ~ Q 0 ) Pk-li • • • 1 Po] 

(3.7) 

Z{t ) = • • • , y, D k ~ l u, . . . , uf 

(3.S) 

We now discuss a way by which the signals u, y, and all necessary 

derivatives of u 


and y in (3.5) may be parameterized using MDFT data. It is well-known that the D FT, 
hence the MDFT, gives an exact parameterization of a signal if the signal is sinusoidal 
and its frequency is exactly equal to one of the DFT frequencies. Specifically, it can be 
shown that, if n is an integer between 0 and N/ 2, 
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(3.9) 


DFT [ j4sin (IS < + ^)] = Mel4> = ^ x p[^-f)] 

at the frequency u n = 2irn/NAt radians/sec, and is zero at all other frequencies. In 
(3.9), A, 9 , M, <j>, N, and At are constants, with N and At defined as in Section 2. 

That is, for the sinusoidal signal defined in (3.9), the magnitude of the time-domain 
signal is given exactly by 2/N times the magnitude of the DFT, and the phase angle of 
the time-domain signal is given exactly by tt/ 2 plus the angle of the DFT. Therefore, if 
the plant inputs are bandlimited and contain sinusoidal components at only the DFT 
frequencies, u(t) and y(t ) can be exactly parameterized in terms of MDFT data. 

For the ideal case assume u(t) is of the form 

N / 2 

“(*) = £ U n sin (u n t + o n ) (3. 10 ) 

where uj n = 2im/NAt and U a n are constants. As stated in Section 2, the frequencies 
W " are the DFT fluencies. Note that, if N is large, u(t) can consist of a large number 
of sinusoidal components. Assuming steady-state has been reached, the plant output 
corresponding to (3.10) is 

N / 2 

y( f ) = £ Y n sin (u n t + (3 n ) (3.11) 

where Y n , /3 n are constants. 

Choose any desired set of k distinct DFT frequencies { W/1 , w /2 , ..., u > Ik ). We call 
these the identification frequencies. Taking MDFTs of u(t) and y{t) at the 
identification frequencies and utilizing the appropriate magnitude and phase 

information can be shown to yield the following parameterizations for u, y, Du, and Du 

[4]: y 



u = C/; 1 sin(a /1 ) + 6 r 72 sin(fv /2 ) -(- . . 

■ +U Ik sin (a lk ) 

(3.12) 


V ^ /i s in(^/|) + F / 2 s in(/?/ 2 ) + . . 

■ + Y ik s ™(Pik) 

(3.13) 

Du 

= ^i\U i l cos{a Il )-\-u} [2 U j 2 cos{a l2 )-\- . 

■ ■ + u >ikU Ik cos(a Ik ) 

(3.14) 

Dy 

= u iv* n cos {0n) + u l2 Y I2 cos(/3 I2 )+ . 

■ ■ +u Ik Y Ik cos(0 Ik ) 

(3.15) 
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where a h are the magnitude and phase of the sinusoidal component in ti (t) at 
identification frequency W/j , and Y h , fi,, are the magnitude and phase of the sinusoidal 
component in y(t) at identification frequency u h . These magnitudes and phases are 
derived from the MDFT as explained above. All higher derivatives of u(f) and y(t) may 
be parameterized in a like manner (see [4]). In this way, all entries in the vectors Y(t) 
and Z(t) in (3.5) can be parameterized with constants at every time step via the MDFT. 
Therefore, the parameters ip can be estimated from (3.5) using standard parameter 
estimation techniques. 

Incidentally, note that, although parameterization for derivatives of u and y are 
obtained, no derivatives of the signals are actually taken. Noise in the signals introduces 
errors in Y{t ) and Z(f), but at worst the errors remain in the same proportion to the 
signal magnitude with higher derivatives and do not increase in size relative to signal 
magnitude, as they would if noisy signals were differentiated. 


Example 

Consider the following 8th-order plant [4]: 



The plant input v(t) is chosen to be a 0.01 Hz, ± 10 V square wave with an 
additive noise component which is uniformly distributed between ± 10 V. For the 
MDFT, we choose N = 5000 and At = 0.02. The identification frequencies are then 
chosen to be 0.01 Hz, 0.03 Hz, and 0.11 Hz, which correspond to DFT frequencies 

and are present in u(t). 

Figure 1 shows parameter estimates for a 2nd-order model found using the 
recursive rational interpolation method described above and applying the recursive 
least squares algorithm to find the parameters from (3.5). The model has the form 


Gr{s) 


CjS ~t~ Cq 

s 2 + a!S + ao 


(3.17) 


where, from Figure 1, Cl = 15.021, c 0 = 4.798, ^ = 5.958, a 0 = 4.795. 

A comparison of G P (s) and G R (s) is shown in Figure 2. There is good 
agreement at all frequencies. It should be noted that u(t) is not bandlimited, as 
specified in (3.10) for ideal inputs, and the noise level in u(t) is comparable to the 
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magnitude of u(t) itself. This indicates that the algorithm is robust to nonideal 
conditions. 


4. ON-LINE MODEL VALIDATION 

It is generally not automatically known whether or not an identified model is 
good, in the sense that its frequency response closely matches that of the plant over a 
range of frequencies. This usually necessitates a model validation stage following the 
identification of a candidate model by any method. 

It has been found in practice that the identification frequencies for interpolation- 
based schemes must be chosen with care in order for the algorithm to find a good 
model. Obviously, choosing correct identification frequencies is difficult if the plant is 
unknown. This indicates the need for some type of model validation before a model 
identified by this method can be used with confidence. 

The identification algorithm presented in the previous section has the built-in 
capability for on-line model validation due to the fact that frequency-domain plant data 
are available via the MDFT. Thus, the algorithm may be given a measure of 
intelligence” by increasing the number of identification frequencies and simultaneously 
deriving several models using different combinations of identification frequencies. The 
frequency response of each model can then be compared with the plant at all 
identification frequencies. The model with the least overall fit error can be taken as the 
best model for the plant of those derived. 


Example 

Let us consider this strategy applied to the example of Section 3, i.e., we 
attempt to find a 2nd-order model for the 8th-order plant (3.16). Since we are 
unsure which identification frequencies to use for the identification, we specify a set 
of frequencies which span the system passband, e.g., choose [ W/1 , W/2 , W/3 , W/4 , W/5 , 
— 27t[ 0.001, 0.01, 0.1, l.o, 10.0, 100.0]. Of course, this necessitates an input 
which is different from the square wave used in the previous example. It can be, 
for instance, a sum of six sinusoids for this example. A model is derived using each 
distinct pair of frequencies, i.e., we derive a model using (u n , u I2 ) as the 
identification frequencies, another model using (u n , o> /3 ), etc. In all, 15 different 
models can be derived taking the above identification frequencies two at a time. It 
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is important to remember that the 15 models are derived simultaneously. 

If we take as an error measurement the distance between the plant and 
model responses in magnitude and phase summed over all identification 
frequencies, we get a figure of merit for each model. Note that these distances are 
easily calculated because plant magnitude and phase data at the identification 
frequencies are available via the MDFT. To illustrate, let us take the following as 

the fit error: 

Error = ^{| 201og 10 (mag(G P (ia; /l )))-201og 1 o(mag(G fi (ja;; l ))) | 

+ | arg(Gp(ju)/,)) - aig{G R {ju>u)) |} (4- 1 ) 

where | ■ | stands for absolute value, mag stands for magnitude, and arg stands for 
argument in degrees. Of course, this definition of fit error is completely arbitrary. 
Whatever the error criterion, however, it must give a single number which 
quantifies the fit error between the plant and model and be capable of being quickly 

calculated. 

The model which yields the lowest error can be taken as the best model of 
those derived. For the present example, we get the following errors correspond,^ 
to the above-defined identification frequencies and error criterion: 


I.D, Freq. Error 
W/ 2 ) lo.585 

(w/i, w/ 3 ) 14.405 

(w/j, w /4 ) 3.6938 
(u> n , w /5 ) 4.2825 
(w/!, uj I6 ) 4.3242 


I.D. Freq. Error 

(u;/2, w/ 3 ) 14.394 

(u>i2, ^'u) 3.6617 
(u) /2 , w /5 ) 4.2390 
(w/ 2 , w /6 ) 4.2805 
(W/ 3 ,W/4) 7 ' 5107 


I.D. Freo. Error 
(u> /3 , ui I5 ) 10.164 
(u>/3, ai/ 6 ) 10.242 
(o>/ 4 , u;/ 5 ) 81.438 
(u>/ 4 , uj / 6 ) 81.370 
(w /5 ,w /6 ) 78.089 


Obviously, the model derived using the frequency pair (w /2 , w J4 ), or ROM #7. 
yields the lowest error and therefore the best fit. This model is given by (3.17) 
with Cl = 17.392, c 0 = 5.315, a a = 7.334, « 0 = 5.320. A comparison of this model with 
the plant is given in Figure 3. Notice that this model is better than the one 
derived in Section 3, since we have “optimized” the identification frequencies via 

on-line model validation. 
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If no model can be found whose fit error is sufficiently small, the model order can 
be easily increased. This amounts to changing the number of entries in Z{t) and 
correspondingly in if. Note that this can be done without changing the inputs or the 

identification frequencies. The order of the calculation (3.5) is simply incremented in 
the algorithm in real time. 

Finally, since the input and output MDFTs are monitored in this method, it is 
possible to make decisions on how fast to update the model parameters or whether to 
update them at all. If the MDFTs are not changing and a suitable model has been 
determined, the parameter estimation part of the algorithm can be shut off. In this 
case, only the MDFTs would proceed as the plant signals evolve. Thus, the monitoring 

of the plant continues, but the parameter estimation part of the algorithm rests if there 
is no need for a parameter update. 


5. AN EXTENSION TO MIMO SYSTEMS 

The above scheme for SISO identification via recursive rational interpolation can 
be extended to MIMO plants as follows. Assume the plant has input vector u(t) = [u v 
« 2 , u,] T and output vector y(t) = [y lt y 2 , yjT The M IMO identification process 
finds constant matrices P, ~ (p x p) and Q { ~ (p x q ), * = 0, 1, m-1 such that the 

model transfer matrix given by 


T(s) = P~\s)Q(s) 

where 

m-l 

P{D) = D m I + jr P t D' 

1=0 

m-l 

Q(D)= £ Q,D' 

1=0 

matches the plant transfer matrix on a given set of frequencies. 

The model can be described by the matrix differential equation 

P{D)y(t) = Q(D)u(t) 
which can be rewritten in the form: 
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D m y - —P m -\D m x y- ■ ■ ■ —PoV + Qm-iD u+ ... +Qo u 

This can be put in a form analogous to (3.5). 

Y(t) = V T Z(t) 


(5.5) 


(5.6) 


where now 


Y{t) = D m y 


$ 


T = [— P m -\i • • • i — Pqi Qm-\1 • ■ ■ 1 ^*>] 


(5.7) 

(5.8) 


and 


Z(t) = [D m_1 y T , • • • , y T , • • • > “1 


(5.9) 


Note that F(t) ~ P x 1, * T ~ P x m{p+q), and Z{t) ~ m(p+y) x 1. 

Again, choose any set of k distinct DFT frequencies {u> n , u; /2 , w /fc }. Denote 
the set of identification frequencies as S = {u n , u I2 , •••, w n-)- Choose integers /, > 0, 
i = l, 5 such that 


i= l 


(5.10) 


Let S i=l, 2, .... 9 be subsets of S such that S, contains I, elements and S.flS^O, 
j J i ' Now the input u,(<) is specified as containing frequency components at the 
frequencies contained in S, The input «, may contain any other frequency componen s 
also, with the exception that the frequency components of S, must be unique to u. 

Again MDFTs are taken of u (t) and y(t) at the appropriate identification 
frequencies. Note that, in the MIMO case, MDFTs are taken of the input u,«) at only 
the identification frequencies contained is S„ but the MDFTs of the output p,(f) must 
be taken at all identification frequencies. 

As in the SISO case, the quantities in Y(t) and Z(t) in (5.6) are parameterized 
from MDFT data at each time instant. Therefore, standard parameter estimation 
techniques may be employed to solve for the parameter matrix * in (5.6). The 
minimum number of identification frequencies necessary for unique identification an 
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the minimum number of identification frequencies per input are topics of current 
research. 

Example 

As a numerical example, the 10-state, 2-input, 2-output plant TGEN given 
in [7] is considered. This is a linearized model of a large turbo-generator. The 
original nonlinear model is given in [8], The A, B, and C matrices for the state- 
space representation of the linearized model are given in [9]. The corresponding 
plant transfer matrix is given by 


G P (s ) 


G Pn (s) G Pl2 (s) 
Gp2i($) G P 22(5) 


(5.11) 


G'p(s) is given in the Appendix. 

A 6-state model with observability indices equal to 3 is chosen for this 
system. Thus, for this example, p = 2, q = 2, m = 3, and k is chosen as 6. The 
identification frequencies are chosen as [w n , u n , w /3 , u l4t w /5 , w /6 ] = 2tt[0.01, 0.1, 
0.2, 0.5, 1.0, 2.0]. The plant inputs are specified as: 

3 

u i(0 = H sina; /^ (5.12a) 

1=1 

6 

*b>(0 = Yl s ' nuj nt (5.12b) 

»=4 

The model transfer matrix identified for this system using rational 
interpolation, G/j(s), is given in the Appendix. Also considered for comparison are 
a 6-state model for this plant given in [9] (denoted as G\f), and a 6-state model 
derived from a balanced realization of the plant (denoted as (7g). These models 
are also given in the Appendix. Figures 4 - 7 show comparisons of the magnitude 
and phase characteristics of G P with those of G R , G M , and G B . 

To get a numerical measure of the closeness of the models to the plant, the 
plant magnitude and phase responses are compared with those of the various 
models at 100 frequencies logarithmically spaced between 10 -3 and 10 1 cycles/sec. 
That is, the following fit error is calculated: 
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(5.13) 


E x = E £ E{l 2010 Slo( ma g( G Pu(>»)))- 201 °Slo( ma g( G ^u<i W n)))l 
i = l j=l n=l 

+ | arg(G P ,j( jw„)) - arg(G XlJ (jw n )) | } 

where X = R, M, and B, and u; n = 2 tt x Kp^"- 1 )/ 99 . This measure is, of 
course, completely arbitrary. However, it is one that is easily and accurately 
calculated and measures ROM fit in the frequency range 10" 3 - 10 1 cycles/sec, 
which is the critical range for this plant. Note that this calculation has nothing to 
do with on-line model validation (Section 4), but is being done simply to compare 
the goodness of the various ROMs. The resulting fit errors are calculated to be 
E r = 3.0157 x 10 3 , E m = 2.4135 x 10 4 , and E B = 9.4244 xlO 3 . Although G R has 
the lowest error, it should be kept in mind that G R is designed specifically to 
match the plant in the frequency range 10 -3 - 10 1 cycles/sec. 


6. CONCLUSIONS 

A method of adaptively identifying reduced-order models for SISO and MIMO 
plants has been presented. The method is based on matching the plant and model 
transfer matrices at a number of frequencies chosen by the designer. The method 
recursively monitors frequency components in the plant inputs and outputs via the 
Moving Discrete Fourier Transform (MDFT). This is a computationally efficient 
method of recursively calculating the DFT of an evolving function of time. The 
identification algorithm has the capability for on-line model validation. 

There is no necessity in this method for the plant input and output signals to be 
purely sinusoidal because the MDFT filters out sinusoidal components in these signals. 
As always, better results are obtained if the system, whatever its order, has a good 
reduced-order model. Computer simulations indicate that the algorithm is robust to 
additive input and output noise and non-bandlimited inputs. 

Areas of current research include error analysis of the algorithm for nonideal 
inputs and stochastic disturbances, and determination of meaningful on-line figures of 
merit for identified models. 
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APPENDIX 


The plant transfer matrix is given by (5.11) with 


<jpn(.s) 


b 7 s 7 + fe 6 .s 6 -f . . . + &o 

A P 


with & 7 = — 4.6075 x 10 1 , b 6 = -4.5705 x 10 3 , b s = - 4.7153 x 10 6 , b 4 - -1.8619 x 10 8 , 
^ 3 = _2.4976 x 10 9 , 6 2 = -1.2507 x lO 10 , = —1-5146 x 10 10 , 6 o = -3.9506 x 10 9 , 


G P i 2 (s) 


b 6 ^ 6 + fr 5 .s 5 + • • • + bp 

A P 


with fe 6 = — 2.0048 x 10 4 , 
6 2 = — 1.2373 x 10 12 , &! = 


b. = -1.7916 x 10 6 , & 4 = - 2.3909 x 10 9 , 


-2.8374 xlO 12 , 6 0 = — 1.8166 x 10 12 , 


. x _ bjS^ + b 6 s 6 + . . . + bp 
^P2i\ s ) — ~aI 


6 3 = — 1.1727 x 10 n , 


with & 7 = -1.1590x10°, 6 6 = —6.9327 x 10 1 , & 5 = 5.6807 x 10 5 , 6 4 = 2.3299 x 10 7 , 
63 = 2.9400 x 10 8 , b 2 = 1 .31 18 x 10 9 , b, = 1 .4696 x 10 9 , b 0 = 1.3869 x 10 9 , and 


Gp 22 {s) 


b s s 8 + bjS 7 + . . . + b 0 
A P 


with 5g=8.1715 x 10 2 , 6 7 = 6.5386 x 10 4 , 6 6 = 7.9770 x 10 7 , fe s = 3.9454 x 10 9 , 

6 4 = 4.6442 x lO 10 , 6 3 = 2.8534 x 10 u , & 2 = 1.8656 x 10 12 , b x = 3.7379 x 10 12 , 

6 0 = 2.0592 x 10 12 . 

In the above, A P = s'° + a 9 s 9 + . . . + a 0 with a 9 = 1.0098 x 10 2 , a 8 = 1.0256 x 10 5 , 
a 7 = 4.2252 x 10 6 , a 6 = 6.5835 x 10 7 , a 5 = 5.4549 x 10 8 , a 4 = 3.2781 x 10 9 , o 3 = 1.4535 x lO 10 , 
a 2 = 2.7412 x lO 10 , a, = 1.8901 x lO 10 , a o = 3.1027 x 10 9 . 

The ROM transfer matrix G R (s) is given by 


Gr n(-s) 


b 5 s 5 + b 4 s 4 -f . . . + b 0 

Ar 


with 6 5 = -7.7158 x lO" 3 , b 4 = 7.6781 x 10~ 2 , b 3 = - 4.6515 x 10\ b 2 = - 4.3011 x 10 2 , 
6 1 = -5.8884 xlO 2 , 6 0 = -1.5861 x 10 2 , 


Gr 22 ( s ) 


b 5 s 5 + b 4 s 4 + . . . + 

Ar 
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with 6 5 = 3.9934 x 10°, 6 4 = -2.6408x 10°, 6 3 = 1.1580 x 10 3 , b 2 = -3.6478 x 10 4 , 
6 X = -9.6222 x 10 4 , 6 0 = -6.4423 xlO 4 , 

G mi (») = *sfl +V 4 + •■■+*« 

n R 

wjth 6 5 = -5.9553 xlO“ 3 , 6 4 = 3.2304 x 10~ 2 , 6 3 = 6.0336 x 10°, 6 2 = 4.6175 x 10 1 , 
6j =5.2247 x 10 1 , 6 0 = 5.5681 x 10 1 , and 


G R22 (s) = M 5 + V 4 + ... + 6 0 

with 6 5 = — 1.64 62 xlO 1 , 6 4 = 1.1952 xlO 3 , 6 3 = 3.1988 x 10 3 , 6 2 = 5.7420 x 10 4 , 

61 = 1.3469 x 10 5 , 6 0 = 7.4649 xlO 4 . 

In the above, A R = s 6 + <z 5 s 5 + ... + a 0 with a 5 = 1.1105 x 10 1 , a 4 = 7.2366 x 10 1 , 
a 3 = 4.5892x 10 2 , a 2 = 1.0098 x 10 3 , a x = 7.4376 x 10 2 , a 0 = 1.2457 x 10 2 . 

The ROM transfer matrix G^v/(<s) is given by 

with 6 5 = 4.8290 xlO" 2 , 6 4 = 1.87 57x 10°, 6 3 = 1.9494 x 10°, 6 2 = -1.2283 x 10 3 , 

6j = -5.76 1 6 x 10 3 , 6 0 = -2.3848 x 10 3 , 

G MI2 ( S ) = V + V + .-.+to 

with 6 S = 2.9153 XlO- 2 , 6 4 = 1.2889 x 10°, 6 3 = -3.0414 x 10' 1 , 6 2 = -2.3880 x 10 4 
b i = -8.2902 xlO 5 , 6 0 = -9.5445 xlO 5 , 

G M21 (s) = b ^±h^±iii±h 

with 6 5 = -4.9194 x 10-\ 6 4 = 9.4826 x 10°, 6 3 = -4.9046 x 10 1 , 6 2 = 5.5 1 64 x 10 2 , 

6j = -9.3964 x 10 2 , 6 0 = 9.0758 x 10 2 , and 


G m (s) -- V 5 + ^ + - 

with 65 = -6.7159 XlO- 1 , 6 4 = 8.0447 xlO 2 , 6 3 = 2.7814 x 10 4 , 6 2 = 8.08 52 x 10 4 , 
6j = 1.2135 x 10 6 , 6 0 = 1.0820 x 10 6 . 
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In the above, A M = s 6 + a 5 s 5 + - - - + <*o with a 5 = 2 - 8238 x 10 1 , a 4 - 2.5831 x 10 2 , 
a 3 = 1 .4679 x 10 3 , a 2 = 8.2145 x 10 3 , a x = 8.8012 x 10 3 , a 0 = 1-6306 x 10 3 . 

The ROM transfer matrix G B {s) is given by 

b^s 6 + b s s s + ■ ■ ■ + bp 

t^BllV 5 / - A g 

with 6 6 = — 1.9102 x 10 -2 , 5, = -6.0759 xlO"\ b 4 = -1.4699 x 10 1 , b 3 = - 4.5056 xlO 1 , 

5 2 = _i.4480 x 10 3 , = -3.5314 x 10 3 , 6 0 = -1.1386 x 10 3 , 

(jBi2f s ; — a b 

with 6 6 =— 1.8202X10-’, h= -4.3626x10-’, b, = -5.6174 x 10”, (,,=3.0908 x 10’, 

(,, = -3.1770 x 10 4 , (,,=—4.7812 x 10 5 , b Q =- 5.2357 x 10 s , 

+ b 5 a 5 + . ■ ■ + fc 0 

with b 6 = 1.5093 x 10- 1 , 6 S = - 2.9259 x 10°, fe 4 = 1.6680 x 10\ 6 3 = -1.1690 x 10 2 , 
b 2 = 6.2777 x 10 2 , = 1.0947 x 10 2 , b 0 = 3.9972 x 10 2 , and 

+ b 5 .5 5 + . . ■ + fr 0 

<JB22t S J - 

with i 6 = 1.6164 xlO" 1 , 6 5 = -9.4448x10°, 6 4 = 1.0597 x 10 3 , f> 3 = 1.5933 x 10 4 , 
i> 2 = 7.2009 x 10 4 , b x = 6.9226 x 10 5 , b 0 = 5.9350 x10 s . 

In the above, A B = s e + a 5 s 5 + • - - + a 0 with a 5 = 2.0310 x 10\ a 4 = 1-6829 x 10 2 , 
a 3 = 9.9166 x 10 2 , a 2 = 4.7178 x 10 3 , a x = 4.8675 x 10 3 , a o = 8.9425xl0 . 
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